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Abstract 

The accumulation point of the period-tripling bifurcation cascade in complex quadratic 
map was discovered by Golberg, Sinai, and Khanin (Russ.Math.Surv. 38:1, 1983, 187), 
and independently by Cvitanovic and Myrheim (Phys.Lett. A94:8, 1983, 329). As we 
argue, in the extended parameter space of smooth maps, not necessary satisfying the 
Cauchy - Riemann equations, the scaling properties associated with the period-tripling 
are governed by two relevant universal complex constants. The first one is 5\ = 4.6002 — 
8.9812z (in accordance with the mentioned works), while the other one, responsible for 
the violation of the analyticity, is found to be 82 — 2.5872 + 1.8067i. It means that in 
the extended parameter space the critical behaviour associated with the period-tripling 
cascade is a phenomenon of codimension four. Scaling properties of the parameter space 
are illustrated by diagrams in special local coordinates. We emphasize a necessity for a 
nonlinear parameter change to observe the parameter space scaling. 

PACS numbers: 05.45.-a, 05.10.Cc, 05.45.Df. 



1 Introduction 



One of the most popular illustrations in the nonlinear science is a picture of the Mandelbrot 
set 0, (Fig. 1). On the complex parameter plane it is defined as the set of points A, at which 
the iterations of the complex quadratic map 

z' = X-z 2 , X,zeC (1) 

launched from the origin z = never diverge. (Here the prime marks the dynamical variable 
relating to the next iteration, i.e. to the next moment of discrete time.) 

The Mandelbrot set has subtle and complicated structure, which is a subject of extensive 
researches. Objects analogous to the Mandelbrot set are present in parameter planes of other 
complex analytic iterative maps too 0, |3|. 

The Mandelbrot set is regarded as a classic example of fractal; this suggests a sort of 
self-similarity, or scaling.^] One particular manifestation of scaling is the Feigenbaum period- 
doubling cascade. The leaves of the Mandelbrot cactus placed along the real axis correspond 
to periods 2, 4, 8,...2 fc ..., and, in asymptotics of large k, they perfectly reproduce each other 
under subsequent magnification by a real universal factor 5p = 4.6692. The limit point of the 
period-doubling accumulation is disposed on the real axis at A^ = 1.401155..., and it is called 
the Feigenbaum critical point. The scaling properties follow from the renormalization group 
(RG) analysis developed by Feigenbaum || Behaviour of the solution of the RG equation 
in complex domain was discussed e.g. by Nauenberg || and by Wells and Overill |J. 

One can find many other critical points in the parameter plane of the map ([!]), at which 
the Mandelbrot cactus displays other scaling regularities. In particular, we can select a path 
on the complex plane across the sequence of leaves corresponding to periods 3, 9, 27..., 3 fc ..., 
and arrive at the period-tripling accumulation point 

A = 0.0236411685377 + 0.7836606508052*. (2) 
1 In fact, as shown by Milnor 0, the subtle structure ("hairiness") of the Mandelbrot set does not reproduce 
itself on deep levels of resolution in small-scales. Speaking more accurately, one should relate the property of 
self-similarity rather to a definite subset called the Mandelbrot cactus || . It consists of a domain of existence 
of a stable fixed point and of domains associated with attractive periodic orbits, which originate from the fixed 
point via all possible bifurcation sequences (see the grey colored part of the picture in Fig. 1). 



As this point was first discovered by Golberg, Sinai, and Khanin [|10|, we call it the GSK critical 
point. Independently, it was found and studied also by Cvitanovic and Myrheim || |TTJ. As 
in the case of period-doubling, the RG analysis has been applied, and self-similarity of the 
sequence of leaves of the Mandelbrot cactus has been revealed. According to Refs. [|K| [TT|, |], 
the scaling factor in the case of period-tripling is complex, 8\ = 4.6002 — 8.9812z.0 

Sometimes it may be preferable to deal not with the one- dimensional complex analytic maps, 
but with the equivalent class of two-dimensional real maps satisfying the Cauchy - Riemann 
equations. 

It would be interesting to discuss a possibility of observation of phenomena associated with 
complex analytic dynamics and the Mandelbrot set in physical systems. Recently this question 
was posed by Beck He has considered motion of a particle in a double- well potential 

in a time-depended magnetic field. As he has stated, under certain assumptions dynamics of 
the particle may be described by the complex quadratic map or by other complex analytic 
maps. It should be emphasized, however, that discussing physical applications for any special 
type of dynamics we should take into account a robustness of a phenomenon under study. In 
this context it is important to know either phenomena associated with the complex analytic 
dynamics survive or not in slightly modified maps, which can violate the Cauchy - Riemann 
conditions. 



In the original work |L0| the authors claimed that the infinite period-tripling cascade could 
occur typically in two-parameter families of real two-dimensional maps (see also [[HJ). In con- 
trast, from studies of Peinke et al. ||H |, Klein ||15|| , and from qualitative analysis of Cvitanovic 
and Myrheim [||] it follows that in non-analytic maps the parameter space arrangement is 
changed drastically, and, apparently, there is no opportunity of surviving for the universality 
intrinsic to the period-tripling. Recently Peckham and Montaldi have presented extensive bi- 
furcation analysis of a complex quadratic map with a non-analytic term [16|, [17| . However, the 
question what happens with the bifurcation cascades of period-tripling remains not clear .[] Sure 

enough, this matter deserves accurate quantitative analysis in terms of the RG approach, and 

2 Beside (|^), a complex conjugate period-tripling accumulation point exists. It is placed at A = A*, and 

scaling in its neighborhood is governed by the conjugate constant <5*. 

3 The same assertion relates to other multi-tupling bifurcation cascades, although we do not touch them in 

the present paper. 



this is the main goal of the present article. Our final conclusion is that in real two-dimensional 
maps the critical behaviour associated with the infinite period-tripling cascade will occur as a 
phenomenon of codimension 4, i.e. it may be found at some special points of a four- dimensional 
parameter space. In other words, the codimension is so high that it seems not easy to observe 
the phenomenon in physical systems, at least in absence of some special symmetries. 

In section 2 we recall contents of the RG analysis developed in fLOl |TT , Then, it is 
extended to study the fixed point solution of the RG equation in respect to a class of per- 
turbations, which can violate the Cauchy - Riemann conditions. We find that among the 
non-analytic eigenmodes of the linearized RG transformation one is of relevance. In section 3 
we formulate the model map appropriate for a study of the extended parameter space. Then 
we explain a necessity of nonlinear variable change to define local coordinates for observation 
of scaling in the parameter space, and construct a procedure to find this coordinate change 
numerically. In section 4 we discuss some details of the parameter space arrangement near the 
GSK critical point and present graphical illustrations of the intrinsic scaling properties. 



2 Renormalization group analysis 

To start, let us follow Refs. ]TU|, [LJ], ||] and assume that discrete time evolution of the complex 
variable z is governed by the complex quadratic map ([!]). Then, for three iterations we have 

z' = f(f(f(z))) = X-(X-(X-zr) 2 . (3) 

Next, we rescale the variable z by a factor «o and choose this factor to normalize the new 
mapping to unity at the origin. Then we have 

= m, (4) 

where 

f 1 (z) = a f(f(f(z/a ))), and a = l//(/(/(0))). (5) 

This is just one step of the RG transformation. 

Now we can apply the same procedure to the function f\(z) and obtain the next function 
/^(-s), and so on. By multiple repetition of the transformation we arrive at the recurrent 



functional equation 

fk+x(z) = a k f k (f k (f k (z/a k ))), a k = 1/ f k (f k (f k (0)))- (6) 

Here f k (z) represents the evolution operator for 3 h iterations of the original map in terms of 
the renormalized dynamical variable. Note that the rescaling constants a k , which appear at 
the subsequent steps of the procedure, are complex. 

According to Refs. |H], [11], f|, at the GSK critical point the sequence f k (z) converges to the 
fixed point function g(z) = lim f k (z), which obeys the equation 

h — >oo 

g(z) = ag(g(g(z/a))) } a = lim a k = l/g(g(g(0))). (7) 

fe^oo 

It is a generalization of the Feigenbaum-Cvitanovic equation || [7| for the case of period-tripling. 
This functional equation may be solved numerically by means of approximating the function 
g(z) via a finite polynomial expansion and computing the coefficients with a use of multi- 
dimensional Newton method [[10], [11]. We have reproduced these calculations, and present the 
results in Table I. Note that only terms of even power are nonzero. 
As follows from the computations, the scaling constant is 

a = -2.09691989 + 2.35827964i (8) 

in good agreement with Refs. JTOj, [H], 

As the rescaling constant is known, it is convenient to redefine the RG transformation. Now 
we will assume that the rescaling factor is the same at all subsequent steps of the procedure 
and equals to the universal value (|[). We just substitute a k = a into Eq.(|6|) and obtain 

f k+ i(z)=af k (f k (f k (z/a))). (9) 

Obviously, this version of the RG transformation possesses the same fixed point g(z), although 
the normalization g(0) = 1 should be regarded now as an arbitrarily accepted additional con- 
dition. 

Let us consider the map z' = g(z) corresponding to the fixed-point solution of the RG 
equation. It may be checked numerically that this map has an unstable fixed point z = z* = 
0.691473 — 0.302692i, and derivative at this point is 

H C = g\z*) = -0.47653180 - 1.054808681 (10) 



Then, as follows from Eq. (0), z*/at is a point belonging to an orbit of period 3. The Floquet 
eigenvalue (or multiplier) for this cycle will be the same complex number \x c . By induction, it 
is easy to prove that there is an infinite countable set of the period-3 fc cycles (z*/a k belongs to 
the period-3 fc orbit), with the same value of the multiplier. 

As the mapping z' = g(z) represents the asymptotic form of the evolution operators for 
the original map ([!]), it must have an infinite set of unstable orbits of periods 3 k at the GSK 
critical point. The same is true for any other map relating to the same universality class. The 
asymptotic value of the multipliers for the set of cycles is the universal constant /i c . 

Now let us assume that we have a real smooth two-dimensional map 

x' = ip(x,y), y' = ip(x,y), (11) 

which depends on some parameters. In terms of the complex variable this map may be ex- 
pressed, in general, via a function of two arguments z = x + iy and z* = x — iy: z' = F(z, z*), 
where F(x + iy,x — iy) = <p(x,y) + iip(x,y). Next, we suppose that at some special values 
of parameters the map becomes complex analytic, i.e. the Cauchy - Riemann equations hold: 
tPx — ipy, <£y — —ipx, and the dynamics is just of the kind intrinsic to the GSK critical point. 
It means that for a number of iterations 3 k (k large) the evolution operator is represented in 
appropriate normalization by the universal function g(z). Now, if we slightly change the values 
of parameters and depart from the original point in the parameter space, some perturbation 
to the fixed point of the RG transformation will appear. In general, the perturbation may not 
satisfy the Cauchy - Riemann conditions, so the evolution operator for 3 fc steps of iterations 
must be written as 

F k (z) = g(z) + eh k (z,z*). (12) 

Here e < 1, and h k is a smooth function of two arguments. By three- fold iteration of the map 
(|l2|) with rescaling of z by factor a we arrive (in the first order in e) to the same form of the 
operator, but with the modified perturbation term: 

h k+1 (z,z*) = a[g'(g(g(z/a)))g'(g(z/a))h k (z/a,(z/a)*) + 
+g'(g(g(z/a)))h k {g{z/a), (g(z/a))*) + h k (g(g(z/a)), (g(g{z/a)))% 

where g' means the derivative of the function g. Obviously, this relation has a structure 

h k+1 = rhh k (14) 



where m is a linear operator defined by the right-hand part of (|I"3|). 

The behaviour of under subsequent iterations of the linear operator is linked with the 
spectrum of eigenvalues. Indeed, let we have an eigenfunction h(z,z*) associated with an 
eigenvalue v. 

uh{z, z*) = a[g'{g{g{z/a)))g'{g{z/a))h{z/a, (z/a)*)+ 

(15) 

+g'(g(g(z/a)))h(g(z/a), (g(z/a))*) + h(g(g(z/a)), {g{g{z J a)))% 
Then, the functional sequence hk may contain a component, which behaves as i> k h(z, z*). Those 
eigenvalues will be of relevance, which are larger than unity in modulus, and which are not 
associated with infinitesimal variable changes. 

To solve the eigenproblem numerically we may expand h(z, z*) over powers of z and z*, and 
use the polynomial representation for g(z), which is already known. As we keep a finite number 
of terms in the expansions, the eigenproblem reduces to a finite-dimensional one. Then, the 
spectrum of eigenvalues may be calculated by means of standard methods of linear algebra. 

As the computations show, the eigenfunction associated with the senior eigenvalue does 
not depend on the second argument z* . The coefficients of the power expansion are listed in 
Table 2; only even-power terms are nonzero. It is clear that this eigenfunction represents a 
perturbation, which retains the map in the class of complex-analytic functions. The computed 
eigenvalue is 

5i = 4.60022558 - 8.98122473z, (16) 

in good coincidence with the data of Refs. [H], [TTj, ^J. 

It appears, however, that one more relevant eigenfunction exists with eigenvalue larger than 
unity in modulus, and the polynomial expansion contains powers of both arguments z and z* 
(Table 3). This solution is responsible for the non-analytic perturbation to the fixed point of 
the RG equation. As we shall see, the eigenvalue may be expressed explicitly via the scaling 
constant a. 



We know that the fixed-point function g(z) is even. Hence, it follows from (15) that 

vh{-z,-z*) = a[g'(g(g(z/a)))g'(g(z/a))h(-z/a,-(z/a)*)+ 
+g > (g{g{z/a)))h(-g{z/al -(g(z/a))*) + h{-g(g{z/a)), -{g{g(z/a)))*)). 
We may represent the eigenfunction as 

h(z,z*) = h s {z,z*) + h a {z,z*), 



;i7) 



where 

h'(z,z*) = h{z > Z * )+ ^- Z >- Z *\ (19) 

is a symmetric part, and 

h*(z,z*) = h{z > Z * ) -^- Z >- Z *\ (20) 
is an antisymmetric component. Then, subtracting fll7l) from ([15]) we have 

uh a (z, z*) = ag'(g(g(z/ a )))g'(g(z/a))h a (z/a, {z/af). (21) 

Now let us set 

h a {z,z*) = ?-^-<$>(z,z*). (22) 

We have g'(z) oc z, so the ratio g'(z)/z does not possess a singularity at the origin. Using the 
relation 

g'(g(g(z/a)))g'(g(z/a))g'(z/a) = 1, (23) 
which follows from ([?]), we obtain a simple equation for Q(z, z*)\ 

v${z, z*) = a 2 $((z/a), (z/a)*). (24) 

Obviously, any product z M (z*) with integer M > 0,N > 0, M + N odd, represents an 
eigenfunction, and the associated eigenvalue is a 2 ~~ M (a*)~ N . Observe that one of these numbers 
is of modulus larger than 1. Indeed, setting M = and N = 1, i.e. $ = z*, we have 

h a (z,z*) = ^-z*, (25) 



and 



Q 2 



v = — = S 2 = 2.58728651 + 1.80679396*. (26) 

a* 



The antisymmetric part of the eigenfunction, as well as the eigenvalue 62 , coincide with those 
obtained numerically (Table 3). 

All eigenvalues beside 5\ and 62 are irrelevant. Some of them are associated with infinitesimal 
variable changes (cf. J7|). They may be expressed explicitly via the fixed point function g (z), 
see Table 4. Numerical calculations show that all other eigenvalues are of moduli less than 1. 

The condition for realization of the critical behaviour in some map is a possibility of tuning 
its parameters to make complex coefficients at two relevant eigenvectors associated with 5± and 



82 equal zero. It yields two complex, or four real, equations on the parameters. Generically, 
the solution will exist if the number of unknowns is not less than the number of equations. So, 
we conclude that the phenomenon may occur as a generic in a space of four real parameters. 
In other words, it is of codimension four. 

In our line of reasoning we have started from the analytic map. Nevertheless, analyticity is 
not a necessary condition for the GSK critical behaviour. To clarify this remark, let us take the 
analytic map g(z) (the fixed point function of the RG equation) and introduce a non-analytic 
perturbation, a small term represented by a smooth function of two arguments, z and z*. It may 
be chosen in such a special way that there will be no contribution from the growing non-analytic 
eigenmode associated with the eigenvalue <5 2 .(See Eq.(p4"|) and the following discussion.) Then, 
under subsequent iterations of the RG transformation the perturbation evolving in accordance 
with Eq. (|T3|) will decay. It means that the evolution operator will tend to the same fixed point 
function as in the analytic map. 

Let us consider a concrete example, the map studied by Gunaratne 13| 



x' = a + x 2 — y 2 , 

(27) 

y' = b + (2 + e)xy. 

By means of variable changes z = x + yi, A = a + ib and e = e/4 this map can be rewritten as 
z' = A + z 2 + e(z 2 — {z*) 2 ). Observe that the non-analytic term (z*) 2 is quadratic, not linear, 
hence, it can not give rise to the eigenvector with eigenvalue 62- As to the contribution into 
the growing analytical eigenmode (eigenvalue Si), it may be compensated by an appropriate 
shift of A = a + ib. Thus, the GSK point must exist; this assertion is in accordance with the 



conclusion of Ref.|T3|- As we have found, for e = 0.2 the critical point is located at a = —0.0355, 
b = 0.7477. 



3 Model map and local scaling coordinates at the GSK 
critical point 

Now we will construct a model map to study a vicinity of the GSK critical point in the extended 
parameter space. 



If we consider a shift of A from the critical value in the complex quadratic map ([I]), then 
only a perturbation given by the eigenfunction h\(z) can arise. As we wish the perturbation 
h,2(z, z*) to appear, it is necessary to add an appropriate non-analytic term. This eigenfunction 
behaves at small \z\ as h 2 oc z*, hence, we rich this goal adding the term proportional to z*. 
Thus, we come to the model map 

z n+1 = \-z 2 n + ez* n , (28) 

which may be regarded as unfolding of the map ([[]). The extended parameter space is the 
two-dimensional complex space C 2 : (A, e). Up to a trivial variable change (z — > —z) this is the 
same map that the authors of Refs. [JT6[, [lTj have chosen for their extensive bifurcation analysis. 
It is equivalent to a four-parameter two-dimensional real map 

x n+1 = a- x 2 n + yl + ux n + vy n , y n+1 = b- 2x n y n + vx n - uy n , (29) 

where we set z = x + iy, X = a + ib, e = u + iv. 

As we have two relevant eigenfunctions hi(z) and h 2 (z,z*), the sequence of the evolution 
operators generated by iterations of the RG transformation will behave as 

f k {z,z*) = g(z) + C 1 (X,e)6 k 1 h^\z)+C 2 (X,e)6 k 2 h^(z,z*y, (30) 

this is true in linear approximation in respect to the deflection from the fixed point solution 
g(z). Here C\ and C 2 are some complex coefficients. They vanish at the GSK critical point 
(A = X c ,e = 0), but become nonzero as we depart from it in the parameter space. 

It would be convenient to use C\ and C 2 as local coordinates in a neighborhood of the 
critical point. In these coordinates the scaling properties of the parameter space are very 
evident. Indeed, let us change G\ to C\jb\ and C 2 to C 2 /5 2 . Then, according to ( |30"D the 
evolution operator fk+i will accept the same form as the operator fk at the old point. It means 
that the similar dynamical regimes will occur at the new point, but with tripled characteristic 
time scale. 

Unfortunately, we do not know explicit expressions of the coefficients C\ and C 2 via the 
parameters of the model map (|28|) . Hence, the coordinate system suitable for demonstration 
of scaling (the scaling coordinates) must be found numerically with a sufficient precision. 



Let us consider a manifold M in the parameter space determined by a condition Ci(A, e) = 0, 
or Ci(a, b; u, v) = 0, where we set A = a + ib and e = u + iv. We can rewrite this equation as 

a — a c + F(u, v), b = b c + G(u,v), (31) 

where a c and b c correspond to the GSK critical point. 

In fact, we may cut the Taylor expansions for F(u,v) and G(u,v) keeping the terms of 
the first and second order. The reason is the following. As we stay at the manifold M, the 
renormalization of the complex coordinate e = u + iv by factor 5 2 has to ensure realization 
of similar dynamical regimes. Let us suppose that we miss some term of power m in the 
expansion of F or C When we perform the scale change e — > s/S^, this term will be of order 
fi-mk ^ g Q ^ err or in the amplitude of the senior eigenvector will behave as 5\/5™ k . In the 
case of our concrete type of criticality (GSK) it is dangerous only for m < 2, otherwise the 
error asymptotically vanishes as k — >• oo. Indeed, in accordance with ( |l6l) and ( p6|) we have 



1 62 1 < and |^2 1 2 < |5i |, but | <5s | 3 > \5i\. So, only terms of the first and second order in the 
Taylor expansion must be evaluated accurately^ and we may write 

F(u, v) = Au + Bv + Pu 2 + Quv + Rv 2 , 

(32) 

G(u, v) = Cu + Dv + Su 2 + Tuv + Uv 2 . 

Here the capital letters A, B,...U designate some constants, which may be found numerically. 

As mentioned, precisely at the GSK critical point the map has an infinite sequence of 
unstable periodic orbits. Let we come out this point and look at some cycle of period N = 3 k . 
In terms of real perturbation vector (x,y) = (Re z, Imz) the evolution over one period is 
determined by the Jacobi matrix 

/ 



^ a ll a 12 



i 0*21 a 22 



(33) 



V 



dx N /dx dx N /dy 
dy N /dx dy N /dy 

While the map is complex analytic, the Cauchy - Riemann equations hold, and an = 022, 
012 = — 021- Otherwise, in presence of the term ez*, these equalities are violated. 
Let us consider trace 

S = a n + a 22 (34) 
4 See other examples of nonlinear parameter change for observation of the multi-parameter scaling in Refs. |2^, 



and determinant 

J = a n a 2 2 — ai 2 a 21 (35) 

of the matrix J at the cycle. Surely, these values depend on the parameters of the map: 
S = S(\,e) = S(a,b,u,v) and J = J(A, e) = J(a,b,u,v). Precisely at the critical point 
(a = a c , h = b c , u = 0, v = 0) they tend, respectively, to 2Re ji c and |/^|, as k grows. 

Let us suppose that we know the first and the second derivatives of the trace and determinant 
in respect to the parameters u =Ree and v =Ime for two sufficiently large periods, 3 fc+1 and 
3 k . We require the following scaling property to hold. If one takes (i) some infinitesimal value 
of e = e = u + iv and (ii) the rescaled value e = S 2 e = (5' 2 u — 5 2 v) + i(S 2 v + S 2 u), where 
5' 2 =Re<5 2 , 5' 2 ' —Im5 2 , then traces and determinants for both cycles must coincide at arbitrarily 
chosen u and v: 

S k +\u,v) = S k (5> 2 u-5>>v,5> 2 v + 5'>u), 

(36) 

J k +\u, v) = J k {5' 2 u - 8%v, 5' 2 v + 5%u). 
Now we substitute the Taylor expansion S(u, v) = S + S u u + S v v + \S uu u 2 + S uv uv + \S vv v 2 , and 
analogous representation for J(u, v), and equalize separately coefficients at the terms of the first 
and of the second power in u and v. It yields ten real equations with ten unknowns A, B, ...U. 
We have found numerically the derivatives (S u , S v , S uu , J u , J v , J uu , ...) for sufficiently large 
periods and solved the mentioned equations. The resulting values of the constants were found 
to be asymptotically stable, estimated as 

A = -0.232750391, B = 0.02699484, C = —B, D = A, 

P = -0.2173, Q = -0.070, R = -0.3184, (37) 
S = 0.1533, T = 0.1010, U = 0.0830. 

(In our computations the maximal periods we could handle were about 3 8 or 3 9 . However, the 
convergence is rather fast, and already at 3 4 3 7 we obtain sufficiently accurate data.) 
Then, the the model map is represented as 

x' = a — x 2 + y 2 + ux + vy + Au + Bv + Pu 2 + Quv + Rv 2 , 

(38) 

y' = b-2xy + vx-uy + Cu + Dv + Su 2 + Tuv + Uv 2 , 

and the coordinate system (a, b, u, v) is appropriate for the demonstration of scaling. Indeed, a 
shift of u and v will contribute exceptionally into the second eigenvector, and a shift of a and 
b — only into the first eigenvector. 



4 Scaling properties of the extended parameter space in 
a neighborhood of the critical point 



The number of real parameters in the model map (pq) is four, and complete bifurcation analysis 
of the parameter space (a, b, u, v) would be complicated and tedious (see, nevertheless, Refs. fll6| , 
I7|). Here we only want to present computer illustrations for scaling properties following from 



the results of Sec. 2. We will consider some two-dimensional cross-sections of the parameter 
space and look at the "topography charts" for dynamical regimes on these surfaces. 

For graphical presentation we will use the technique of Lyapunov charts (see |2D, |2T| for 



previous applications of this method). As we have chosen the parameter-space cross-section to 
be studied, we compute the senior Lyapunov exponent for the model map at each pixel of the 
two-dimensional plot, and mark this pixel by an appropriate grey tone. We code the Lyapunov 
exponent from large negative values to zero by tones from dark to light grey (it corresponds 
to periodic regimes). White color represents zero Lyapunov exponent (e.g. quasiperiodicity or 
states at the onset of chaos), and black — positive values (chaos). Such a convention ensures 
a clear vision of the border between regular and chaotic dynamics. In some domains of the 
parameter space the model map manifests divergence; these areas are colored uniformly by one 
special grey tone. 

In our computations the trajectories are launched from the origin z = 0, but the calculations 
for evaluation of the Lyapunov exponent start after some sufficient number of iterations to 
exclude the transients. 



To begin, let us illustrate the kind of scaling that follows from results of Refs. [[10] , [II], |J. 
In the top part of Fig. 2 (a) we show a cross-section of the parameter space by surface u = 0, 
v = 0, at which the equation Q3"B] ) reduces to a conventional complex analytic quadratic map 
(|H). What is observed is just a fragment of the Mandelbrot cactus. The GSK critical point 
is located exactly at the center of the diagram. We select a small box containing the critical 
point, enlarge and rotate it in accordance with multiplication by 5i = 4.6002 — 8.9812i, and 
then reproduce on the bottom plot. In a course of this procedure we redefine the legend for 
the Lyapunov exponent: as the characteristic time scale triplicates, we decrease by 3 times all 
the values separating the intervals for coding by definite grey tones. Observe excellent visual 



similarity of both pictures. 

Next, let us take another cross- section of the parameter space of the map (|38|) , by the surface 
(a = 0,6 = 0). It means that in the parameter space of the original map (p8|) we are at the 
manifold M: C*i(A, e) = 0. The Lyapunov chart is shown in Fig.2(b); the critical point is placed 
precisely in the middle of the plot. Observe that visually the picture has nothing common with 
the familiar Mandelbrot set. However, the self-similarity does hold, although it is governed now 
by the new scaling constant (^). To demonstrate it, we take a box around the critical point. 
After magnification and rotation corresponding to multiplication by 82 — 2.5872 + 1.80672, and 
with redefinition of the grey-scale coding rule for the Lyapunov exponent, we obtain the chart 
shown on the bottom plot. It looks remarkably similar to the top diagram. 

To give clearer impression of the parameter space arrangement on the manifold M we present 
in Fig. 3 a schema explaining some details of the bifurcation structure. The largest area in the 
middle corresponds to a stable period- 1 state. Here the map possesses the fixed point with 
two complex conjugate multipliers of modulus less than unity. From the bottom this domain 
is bounded by a curve, at which the Neimark bifurcation occurs. Here two complex multipliers 
cross the unit circle. Along the border argument of complex multiplier is varied. For irrational 
arguments, measured in 2tt units, an attractor, which appears as a result of the bifurcation, 
will be a torus. For rational values it is a resonance cycle (at the sharp ends of the Arnold 
tongues). At the upper border of the stability domain period-doubling bifurcation takes place. 
It may be followed either by secondary period doubling (top left part of the border for the 
period-2 domain), or by Neimark bifurcation (top right part of the border). Note overlap of 
the period-1 and period-2 areas, which indicates presence of multistability (coexistence of two 
attractors with distinct basins). At this region the borders of the stability domains appear to 
be associated with the hard bifurcations (jumps). At the left bottom part of this structure 
another similarly arranged smaller formation is attached. Here the period-3 regime occurs in 
the middle part of the stability domain. According to results of the RG analysis and scaling 
arguments, the sequence of the smaller domains, each of which is attached to its precursor, 
must be infinite. Asymptotically they become self-similar, accepting a universal form. (In fact, 
already on the very crude resolution scale of Fig. 3 the domains of periods one, three, and nine 
look remarkably similar.) 



The present illustrations support our basic statement that the second eigenvalue 62 is of 
relevance for the non-analytic modification of the map. The scaling properties of the extended 
parameter space are determined by two constants, the earlier known Si, and the novel 62 ■ 

5 Conclusion 

In this article we discuss scaling behaviour of complex maps at the accumulation point of 
period-tripling in the extended parameter space. 

We consider the fixed point function of the RG equation of Golberg, Sinai, and Khanin, 
and analyze a class of perturbations, which are supposed to be small and smooth, but not 
necessarily obey the Cauchy-Riemann conditions. There is no need in modification of the 
fixed point equation, but we revise the linearized RG equation, which describes behaviour in 
a vicinity of the analytic fixed point. Then, we study the eigenvalue spectrum and detect two 
relevant eigenvectors, one of which is analytic, and another one is non-analytic and violates the 
Cauchy - Riemann conditions. 

To observe the period-tripling universal scaling behaviour in a non-analytic map one needs 
to make two complex coefficients at the relevant eigenvectors equal zero by means of appropriate 
selection of parameters. Typically, it requires controlling not less than four real parameters. It 
may be thought that the possibility of observation for such type of behaviour in physical systems 
is rather problematic. Nevertheless, we can not exclude that the period-tripling criticality may 
be found in some systems, which possess special symmetries. Apparently, this is an interesting 
direction of search for physical applications of the complex analytic dynamics. 
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Table 1: The coefficients of polynomial expansion for the fixed-point function g(z) 



1 


1.0 + 0.0i 


z 2 


0.054665304 - 0.749020944i 


z 4 


-0.024397241 - 0.052466461i 


z 6 


-0.002529112 - 0.001197430i 


z 8 


-0.000088081 + 0.000137556i 


z io 


0.000000729 + 0.000018289i 


z 12 


0.000000541 + 0.000001194i 


z 14 


0.000000074 + 0.000000048i 



Table 2: The coefficients of polynomial expansion for the senior eigenfunction h^\z) 



1 


1.0 + 0.0i 


z 2 


0.181223377 - 0.141361034i 


z 4 


0.002303322 - 0.015962515i 


z 6 


-0.001131785 - 0.000491487i 


z 8 


-0.000151737 + 0.000053131i 


z io 


-0.000010514 + 0.000009686i 


z 12 


-0.000000310 + 0.000001085i 


z 14 


0.000000033 + O.OOOOOOlOOi 



Table 3: The coefficients of polynomial expansion for the eigenfunction h^ 2 \z,z*) associated 
with the eigenvalue 6 2 = 2.5872865f + f .80679396? 





1 


z* 


(z*) 2 








f 


1.0 + 0.0i 


-3.029846 


0.398114 


-0.022670 


-0.002329 


-0.000073 






+6.351031* 


+0.008514* 


-0.073276i 


-0.003512i 


-0.0002278i 


z 2 


f.094504 


0.068249 


0.086111 


-0.001995 


-0.000053 


0.0000004 




-1.043125? 


+1.082134* 


+0.039504i 


+0.001278i 


+0.000130i 


+0.000006i 


z 4 


-0.048082 


0.053978 


0.001900 


-0.000124 


0.0000003 






-0.141290i 


+0.057211i 


-0.001493i 


+0.000212* 


+0.000009i 




z e 


-0.011263 


0.004950 


-0.000450 


-0.000005 








-0.005046i 


-0.003601i 


-0.000460i 


+0.000022i 






z 8 


-0.000616 


0.000280 


-0.000058 


0.0000002 








+0.000710i 


-0.000811i 


-0.000028* 


+0.000002i 








-0.000002 
+0.000109i 


-0.000004 
-0.000074i 


-0.000004 
+0.0000002i 









Table 4: Some eigenfunctions of the eigenproblem (|15|) associated with the infinitesimal variable 
changes 



Variable change, 
A < 1 


Eigenfunction 

h 2 {z,z*) 


Eigenvalue v 


Shift: 

z -> z + A 




v = a = -2.0969 + 2.3583z 


Scale change 
z -> z(l + A) 


zg\z) - g(z) 


v = 1 


Non-analytic change 
z -> z + z* ■ A 


z*g'(z) - (g(z)Y 


v = a/a* = -0.1169 - 0.9931i 
\u\ = 1 




Figure 1: Mandelbrot set for complex analytic quadratic map ([l]). A unification of grey colored 
areas is the Mandelbrot cactus. Periods associated with the leaves of the cactus are marked 
by respective numbers. Two selected boxes show domains of the Feigenbaum period-doubling 
scaling and of the period-tripling scaling. 




to ft) 

Figure 2: Cross-section of the parameter space by a surface u = 0, v = 0, at which the map 
(|38|) is complex analytic (a), and by the manifold M, at which the contribution of the senior 
eigenvector into a perturbation of the RG fixed point is excluded (b). Grey tones from dark to 
light code values of the Lyapunov exponent from large negative to zero. White color corresponds 
to zero, and black - to positive Lyapunov exponent. Divergence is shown by uniform coloring 
with one special grey tone. The GSK critical point is located exactly at the center of the 
diagrams. A small box containing the critical point, is magnified and rotated in accordance 
with multiplication by 8\ ~ 4.6002 — 8.9812z on the panel (a) and by 82 — 2.5873 + 1.8068? on 
the panel (b). The coding rule on the bottom plots is redefined to reveal the self-similarity: all 
the values separating intervals for coding by definite grey tones are decreased by 3 times. 




Figure 3: A schema of the parameter space arrangement in the cross-section of the extended 
parameter space by the manifold M, at which the contribution of the senior eigenvector into 
a perturbation of the RG fixed point is excluded. Bifurcation lines are marked as pd (period- 
doubling), t (tangent), N (Neimark). Light grey areas correspond to periodic regimes of periods 
marked by the Roman figures. Dark grey designates the area of quasi-periodicity and Arnold 
tongues (not shown). On the Neimark bifurcation curve some points of rational arguments of 
the multiplier are marked, these are the sharp ends of the respective Arnold tongues. The GSK 
critical point is marked by a cross. 



